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Abstract 


Recent determinations of the equation of state of dark energy hint that this may well be of the 
phantom type, i.e., Wde < — 1- If confirmed by future experiments, this would strongly point to the 
existence of fields that violate the dominant energy condition, which are known to present serious 
theoretical difficulties. This paper presents an alternative to this possibility, namely, that the 
measured equation of state, Wde-, is in reality an effective one, the equation of state of the quantum 
vacuum, wx = — 1, plus the negative equation of state, Wc, associated to the production of particles 
by the gravitational field acting on the vacuum. To illustrate this, three phenomenological models 
are proposed and constrained with recent observational data. 
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I. INTRODUCTION 


Recently, the analysis of a well of data provided by the Planck satellite [l| strengthened still 
further our conhdence in the, so-called. Lambda cold dark matter (ACDM) model. Thus 
far, it constitutes the most promising cosmological model in the market because, notwith¬ 
standing its simplicity, it hts rather well most observational data. By assuming a spatially 
flat, homogeneous, and isotropic universe, of which the main sources of gravity at present 
are pressureless matter (baryonic plus dark) and the energy of the quantum vacuum (the 
vacuum pressure is related to the latter by p\ = —Pa), h successfully describes, with just six 
free parameters, the evolution of our Universe up to the present era of accelerated expansion. 


A crucial quantity in models based on Einstein gravity aimed at accounting for the current 
accelerated phase of expansion is the equation of state (EoS) of dark energy (DE), the 
ratio Wde = Pde/Pde of fh® pressure to the density of dark energy. The latter drives the 
acceleration thanks to its highly negative pressure, pde < —pde/^- In the case of the ACDM 
model, this agent is nothing but the energy of the vacuum whence the corresponding EoS 
parameter is just w\ = —1. In spite of the success of this model, recent model-independent 
measurements of Wde seem to favor a slightly more negative EoS (see, e.g.. Refs. 3 , 3 , 3 , 
and 3 ), which, if conhrmed, would invalidate the model. In particular, the Planck mission 
yields Wde = 3 - 


Rest et ai, using supernovae type la (SN la) data from the Pan-STARRSl Medium Deep 
Survey in the redshift interval 0.03 < 2 ; < 0.65 found Wde = “l-166lo;ogg, "^de < —1 at 
2.3(7 conhdence level 3- However, they caution that it is unclear whether the tension with 
the ACDM value arises from new physics or a conjunction of chance and hidden systematics. 
Similarly, Cheng ef al. obtained Wde = —1.16 ± 0.06 {wde < —1 at 2.6a) 3- Likewise, the 
authors of Ref. (5|, using geometrical data from SN la, baryon acoustic oscillations (BAOs), 
and the cosmic microwave background (CMB) radiation, determined —1 at ~ 1.9a 

level. They conclude that, at 2a, either the Supernova Legacy Survey [61] and Panoramic 
Telescope and Rapid Response System data 3 contain unknown systematics or the Hubble 
constant is lower than 71 km/s/Mpc, or else, Wde < w\. 
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The simplest way to achieve Wde < —1 in a manner consistent with general relativity is to 
assume that the dark energy corresponds to some or other scalar held, called a “phantom 
held,” that violates the dominant energy condition by allowing its kinetic term to have the 
“wrong” sign js]. Thus, its energy density and pressure are given by pph = —+ ld(0) 
and Pph = — I 0^ — respectively; thereby, its EoS obeys Wph = Pph/Pph < — 1- 


and 


While this kind of helds shows compatibility with observational data -see, e.g.. Ref. |l|- 
constitutes a serious contender of the ACDM model, its blatant violation of the said energy 
condition entails serious problems on the theoretical side. Since their energy density is not 
Dounded from below these helds suher from instabilities at the classical and quantum levels 


id] and other maladies jlll-ll3| that cast doubts about their very existence. Nevertheless, 


observationally, they cannot be discarded. Then, it is natural to wonder whether an ehective 
EoS less than —1 can be accomplished without resorting to phantom helds. The target of 
this paper is to explore whether the combined ehect of the pressure of the quantum vacuum. 
Pa = —PA, and the negative creation pressure of particles from the gravitational held of 
the expanding Universe acting on the vacuum can achieve this. In this new scenario, the 
ehective EoS comes to be Weff = w\ + Wc < —1, where Wc is the EoS related to the creation 
pressure. If the answer is in the affirmative, then the need to recourse to phantom helds 
will weaken. 


This paper is organized as follows. The next section briehy sums up the phenomenological 
basis of particle creation in expanding homogeneous and isotropic, spatially hat, universes. 
Section HTTI proposes three diherent phenomenological expressions for the particle production 
rate. Section |IV] specihes the various sets of data and the statistical analysis employed to 
constrain the models. The corresponding results are presented in Sec. |Vl Section |VT] briehy 
explores whether coupled models of DE may give rise to ehective EoS of phantom type. The 
concluding section summarizes and gives comments on our hndings. As usual, a subindex 
zero attached to any quantity means that it must be evaluated at present time. 
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II. COSMOLOGICAL MODELS WITH PARTICLE CREATION 


As investigated by Parker and collaborators the material content of the Universe may 
have had its origin in the continuous creation of radiation and matter from the gravitational 
held of the expanding cosmos acting on the quantum vacuum, regardless of the relativistic 
theory of gravity assumed. In this picture, the produced particles draw their mass, momen¬ 
tum and energy from the time-evolving gravitational background which acts as a “pump” 
converting curvature into particles. 


Prigogine 


15| studied how to insert the creation of matter consistently in Einstein’s held 


equations. This was achieved by introducing in the usual balance equation for the number 
density of particles, {nu°').a = 0, a source term on the right-hand side to account for 
production of particles, namely. 


= nP, 


( 1 ) 


where is the matter huid 4-velocity normalized so that m" Mq, = 1 and P denotes the parti¬ 
cle production rate. The latter quantity essentially vanishes in the radiation-dominated era 
(not to be considered in this paper) since, according to Parker’s theorem, the production of 
particles is strongly suppressed in that era |^. The above equation, when combined with 
the second law of thermodynamics naturally leads to the appearance of a negative pressure 
directly associated to the rate P, the creation pressure Pc, which adds to the other pressures 
(i.e., of radiation, baryons, dark matter, and vacuum pressure) in the total stress-energy 
tensor. These results were subsequently discussed and generalized in Refs, [l^, [l^, and 


2d, 


191 by means of a covariant formalism and further conhrmed using relativistic kinetic theory 


2l|. 


Since the entropy flux vector of matter, ncrw", where cr denotes the entropy per particle, 
must fulhll the second law of thermodynamics {nau°').^a ^ 0, the constraint P > 0 readily 
follows. 


For a homogeneous and isotropic universe, with scale factor a, in which there is an adiabatic 
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process of particle production from the quantum vacuum, it is easily found that [l7|, llQl ] 


P + P-r 


( 2 ) 


Therefore, being pc negative it can help drive the era of accelerated cosmic expansion we 
are witnessing today. Here p and p denote the energy density and pressure, respectively, of 
the corresponding fluid; H = a/a is the Hubble factor; and, as usual, an overdot denotes 
differentiation with respect to cosmic time. Since the production of ordinary particles is 
much limited by the tight constraints imposed by local gravity measurements 22l-l24l| . and 
radiation has a negligible impact on the recent cosmic dynamics, for the sake of simplicity, 
we will assume that the produced particles are just dark matter particles. 


III. COSMIC SCENARIO 

Let us consider a spatially flat Friedmann-Robertson-Walker universe dominated by pres¬ 
sureless matter, baryonic plus dark matter (DM), and the energy of the quantum vacuum 
(the latter with EoS p\ = —Pa) in which a process of DM creation from the gravitational 
held, governed by 

Pdm T Pd^ pdm L (3) 

is taking place. In writing the last equation, we used Eq. ([1]) specialized to DM particles 
and the fact that pdm = ^dm ni, where m stands for the rest mass of a typical DM particle. 
Since baryons are neither created nor destroyed, their corresponding energy density obeys 
Ph + 3Hpf, = 0. On their part, the energy of the vacuum does not vary with expansion, 
hence pa = constant. 


In this scenario, the total pressure is pa + Pc, and thereby the effective EoS is just the sum 
of the EoS of vacuum plus that due to the creation pressure, 

r 


_ Pa , Pc _ , 

Weff - — + — - -1 

Pa Pdm 


( 4 ) 


Since, by the second law, T is positive semidehnite, we have that the effective EoS can 
be less than —1 without the need of invoking scalar fields with the wrong sign for the 
kinetic term. Therefore, thanks to the combined effect of the vacuum and creation pres¬ 
sures one can hope to get rid of phantom fields and of the severe drawbacks inherent to them. 
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The Friedmann equation for this scenario is 


— 2 {Pb + Pdm + Pa)- (5) 

To go ahead, an expression for the rate T is needed. However, the latter cannot be ascer¬ 
tained before the nature of dark matter particles be discovered. Thus, in the meantime, we 
must make ourselves content with phenomenological expressions of T. Here, on grounds of 
simplicity, we assume three phenomenological Ansatze, namely. 


F = 3/3 Id 

(Model I), 

( 6 ) 

F = 3/3 FF [5 — 5 tanh(10 — 12a)] 

(Model H), 

( 7 ) 


and 


T = 3/3 FT [5 — 5 tanh(12a — 10)] (Model HI), (8) 

where /3 is a constant parameter satisfying 0 < /?. As indicated, models I, H, and HI 
correspond to Ansatze ([6]), ([Tj), and ([8|), respectively. 


The ratio T/SiF is a constant in model I. Figure [H shows the said ratio in terms of the scale 
factor for models H and HI assuming /3 = 0.1. It should be noted that for /3 larger than 0.1 
one is led to tCe// < ~2, which lies much away from the reported Wde values. In all three 
cases, F/3iF < 1 at any scale factor. 

Inserting Eqs. ([6]), ([7]), ( 0 ), in ([3]), and integrating, we have 

pdm = Pdmo (9) 


for model I and 


-3 


pdm = pdmo a exp 



P(a) 


da 


( 10 ) 


for the two other models, with g{a) =5 — 5 tanh(10 — 12a) and g{a) = 5 — 5 tanh(12a — 10) 
for models H and HI, respectively. In both cases, g{a)/a —?• 0 as a —)• cx), whence the integral 
stays hnite in that limit. 


In terms of the redshift, z = a ^ — 1, the Hubble expansion rate reads 


— — ^bO (1 + ^)^ + ^dmO (1 + ^)^*'^ + ^AO 


( 11 ) 
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Figure 1. Evolution of ratio r/3-ff as a function of the scale factor for the model II (solid, black 
line) and model III (dashed, red line). In drawing the graphs we have taken f3 = 0.1. 

for the model I and 

— J^2 ^ ~ (1 + + ^dmO (1 + z)^ exp ( — 3/3 f ——+ riAO) (12) 

Ho V Jo {1 + z) J 

for models II and III, where the hljo denote the cnrrent fractional densities of baryons, dark 
matter, and vacnnm, and we have taken into acconnt that the baryon energy density scales 
as (1 + z)^. Here and thronghont the scale factor is normalized so that Oq = 1. 


IV. OBSERVATIONAL CONSTRAINTS 


To constrain the free parameters 9i = {(3, fldmo) of the models above, we nse 580 Snpernova 
type la data points in the redshift interval 0.015 < ^ < 1.41 of the Union 2.1 compilation 


25l |. 9 gamma-ray bnrsts in the redshift interval 1.54 < z < 3.5 j27|, 6 points of baryon 


aconstic oscillations in the redshift interval 0.106 < 2 ; < 0.730 [28j, and 28 data points of the 


Hnbble rate in the redshift interval 0.09 < z < 1.75 


29|. Their best fit valnes, with their 


corresponding la nncertainties, are presented in SEct. El These follow from minimizing the 
likelihood fnnction L oc exp(-xLaz/2) with xLaZ = x|zv/a + Xgrb + xIao + where each 
Xi (specified below) qnantifies the discrepancy between theory and observation. We not nse 
the popnlar CMB shift parameter dne to its comparative heavy dependence on the standard 
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cold matter and ACDM models. 


A. Supernovae type la 

Data from SN la are an important tool for understanding the recent evolution of the Universe. 
Here, we use the Union 2.1 compilation 2^, available at http://supernova.lbl.gov/Union, 
that contains 580 SN la data in the redshift range 0.015 < ^ < 1.41. The distance modulus 
predicted for a given supernova of redshift ^ can be expressed as 


IJ-th{z) = 5 log 


10 


\10pc 


+ ho , 


( 13 ) 


where = (1 + 2 ;) dz' is the Hubble-free luminosity distance and /io = 42.384 — 
51ogioh, with h = iho/(100Km/s/Mpc), the reduced Hubble constant. 


Assuming the SN la data follow a Gaussian distribution, we have 


,2 (n N ^iiho)]^ 

XsNIaK^h hoj — 2_^ 
i=l 




(14) 


where denote the observed value, denote the value predicted by the model, = 
(/3, ^dmo), and cTj stands for the la uncertainty associated to the Ah data point. To eliminate 
the effect of the nuisance parameter juq, which is independent of the data points and the 
data set, we minimize the right-hand side of the last equation with respect to juo following 
the procedure of Ref. (SOj. We hrst expand in terms of /xq as 


— A — fiQ B -j- fiQ C, 


where 


580 


m) = E 


dj, /Up = 0 )]^ 


B[e,) = f] 


i=l 


2 = 1 


and 


580 




1 


2=1 ^ ' 

Equation (fT5|) presents a minimum for no = B/C at 

B\ei) 


x\e,) = A{e,) 


c 


( 16 ) 

( 16 ) 

( 17 ) 

( 18 ) 

( 19 ) 
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Thus, rather than minimizing we minimize x^{Oi) that is independent of po- 

Clearly, Xmm = Xmin- Thus the Hubble constant, Hq, does not enter the calculation of the 
cosmological parameters when using SN la data. 


In our analysis we have not taken into account the small correlations between the SN la data 


points. As seen in the paper by Conley et al. 


the correlations induced by systematic 


errors have a very minor impact (see, e.g.. Fig. 11 in that paper). This was also found by 


Ruiz et al., 


3l| (see, e.g.. Fig. 11 there). We feel therefore conhdent that they will not 


signihcantly alter the overall result of our analysis (i.e., that, as argued below, models I and 
II may well explain the value of less than -1 reported in Refs. [2-5] for the EoS of dark 
energy.) 


B. Gamma-ray bursts 


Gamma-ray bursts (GRBs) are very energetic astrophysical outbursts usually at higher 
redshifts than SN la events. Unfortunately, very often, GRBs are not standard candles as 
is the case of SN la. Therefore, it becomes necessary to calibrate them if they are to be 
employed as useful distance indicators. The uncertainties in the observable quantities of 
GRBs are much higher than in SNs la, since so far there is not a good understanding of their 
source mechanism. These issues favor a controversy over the use of GRBs for cosmological 
applications [26 1. 


Recently, in Ref. 


27|, a set of 9 long gamma-ray bursts (LGRBs) in the redshift interval 


1.54 < .^ < 3.5 was calibrated through the type I fundamental plane. This is dehned by the 
correlation between the spectral peak energy Ep, the peak luminosity Lp, and the luminosity 
time Tl = Eiso/Lp, where Eigo is the isotropic energy. This calibration is one of the different 
proposals to calibrate GRBs in an model cosmological-independent way. The fact that a 
control of systematic errors has been carried out to calibrate these 9 LGRBs 27(| makes this 
compilation compelling. 


9 








The function for the GRB data reads 


XcRBi^i^ /^o) — 


i=l 




( 20 ) 


Here, we perform the same procedure described in the previous section. 


C. Baryon acoustic oscillations 

These can be traced to pressure waves at the recombination epoch generated by cosmological 
perturbations in the primeval baryon-photon plasma and appear as distinct peaks in the 
large-scale correlation function. 


For the BAO measurements, we use the six estimates of the BAO parameter 

r^(z) 1 1/3 


A(z) — ■\/ -j- 


z^F(z) 


( 21 ) 


given in Table 3 of Ref. 2^, that are in the redshift range 0.106 < ^ < 0.730, which is the 
Ho independent. In this expression, r is the comoving distance. 


The function for the BAO data is 


i=l 




( 22 ) 


Again, we neglect the correlations in the BAO data as their impact in the hnal results is 
expected to be rather small (see, e.g.. Fig 11 in Ref. |3l|). 


D. History of the Hubble parameter 


The differential evolution of early type passive galaxies provides direct information about 
the Hubble parameter, H{z). An updated compilation of 28 data points H{z) lying in the 
redshift interval 0.09 < < 1.75 can be found in Ref. [29|. Here, we use these data to 

constrain the cosmological free parameters 9i of the three models under consideration. 
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We compute the Xh function defined as 


xUe., ff„) = f' 

i=l 


cU^i) 


( 23 ) 


where Hth{zi, HQ,9i) is the model-predicted value of the Hubble parameter at the redshift 
Zi- This equation can be recast as 


xU»i, ffo) = Hi f; - 2H„ f; 


<t? 


(24) 

i.i °i s 

The Xh function depends on the model parameters 6i as well as on the nuisance parameter 
Ho, the value of which is not very well known. To marginalize over the latter, we assume that 
the distribution of Hq is Gaussian with standard deviation width hho and mean Hq- Then, 
we build the posterior likelihood function Lh{0) that depends just on the free parameters 




where 




7^h{Ho) = 




(25) 


(26) 


y^/^TTCT Hq 

is a prior probability function widely used in the literature. For Hq, we take the best-fit 


value provided by Riess et al. 


32| . Finally, we minimize Xnidi) = —21nLj:/(6'j) with respect 


to the free parameters Oi to obtain the best-fit parameter values. 


V. RESULTS 

Tables U and mi summarize the main results of the statistical analysis carried out using 
the set of data, SNIa-fGRB-l-BAO and SNIa-|-GRB-|-BAO-|-Ff( 2 ;), respectively. 

Figure [2] shows the 68% and 95% conhdence contours in the plane (3 — VLdmo for model I 
considering the sets of data SNIa-fGRB-l-BAO (left panel), and SNIa+GRB-|-BAO-|-if( 2 :) 
(right panel). Although the best £t for (3 is small for the two cases, the possibility of a small 
creation rate (/9 > 0 ^ F > 0) is not ruled out. In fact, 0 < /9 < 0.093 and 0 < /3 < 0.136 
at la and 2a confidence levels, respectively, for SNIa-|-GRB-|-BAO, and 0.039 < (3 < 0.108 
(0.017 < /3 < 0.130) in la (2a) when we considered SNIa-|-GRB-|-BAO-|-i7(z). We see that 
a nonzero creation rate (F > 0) is consistent with the set data SNIa-fGRB-|-BAO-|-i7(z) 
in la and 2a, and from this analysis, we obtained Weff{z = 0) = —l.OTOlgiQgg in la of 
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statistical confidence. For the joint analysis SNIa+GRB+BAO, we have an effective EoS, 
to,„(z = 0) = -1.02491S:S. This hints that the combined effect of the quantum vacuum 
and the particle creation rate may indeed result in an effective EoS lower than —1. The 
associated uncertainty of tCg// was determined by using the standard error propagation 
method. 



P 

RdmO 

xLJdof 

Model I, Eq. Q, 

U.UZO_Q Qg5 

n 907+0-016 

u.zo/_0 015 

0.9589 

Model II, Eq. ([71), 

0 n‘^+0-081 

D 2'il +0-016 

U.ZOi_Q QJ^g 

0.9587 

Model III, Eq. ([8]), 

U.UUO_Q QQg 

f) 2^11 +0-017 
U.ZOi_o 016 

0.9590 


Table I. Best-fit values with their la errors of the free parameters, /3 and ^dmo, for the three models 
considered in this work obtained from joint analysis SNIa-|-GRB-|-BAO. 



P 

RdmO 

xLin/dof 

Model I, Eq. ([6]), 

0.0729l‘]:°il 

n 9Q9+0.016 
W.ZOZ_o 015 

0.9549 

Model II, Eq. (0, 

0.0158lO:“S2 

n 9QO+0.015 
U.ZOO_o_oi4 

0.9774 

Model III, Eq. ([8]), 

n ni 97+0.0049 
U.U1Z^_0 0050 

0 

0.9705 


Table II. Best-fit values with their Icr errors of the free parameters, (3 and UdmO^ for the three 
models considered in this work obtained from the joint analysis SNIa-|-GRB-|-BAO-|-R(z). 

Figures [3] and m show the 68% and 95% conhdence contours in the (3 — fldmo pfon for models 
II (left panel), and the evolution of Wg// in la (right panel) for the respective data sets, 
SNIa-fGRB-l-BAO, and SNIa-fGRB-l-BAO-l-if(z). As can be observed, the inclusion of the 
data set, H{z), signihcantly reduces the constraints on the parameter (3, since it, changed 
from 0 < (3 < 0.111 to 0.0081 < (3 < 0.024 in la, for example. This scenario, presents as 
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Figure 2. Left panel: la and 2a confidence contours for model I obtained from joint analysis 
SNIa+GRB+BAO. Right panel: The same but using the set of data SNIa+GRB+BAO+Lf( 2 ;). 

effective EoS, Weff{z = 0) = —1 . 2944 lg;gQ 49 for the analysis with SNIa+GRB+BAO, and 
Weff{z = 0) = —1.555lQ;Qgg for SNIa+GRB+BAO+iL(z), both in la. This model, which is 
characterized by a prodnction rate of particles that grow thronghont cosmic history, today 
presents a signihcant effective EoS phantom. Figures [5] and [6] show the 68% and 95% conh- 
dence contours in the /3 — fldmo pla-n for models III (left panel) and the evolution of Weff in la 
(right panel) for the respective data sets, SNIa+GRB+BAO and SNIa+GRB+BAO+iL( 2 :). 
Here, we have that the model presents a small effective EoS in the present moment, being 
Weff{z = 0) = —llo.oois SNIa+GRB+BAO and 

Weffiz = 0) = -1.0021°;°°^ for SNIa+GRB+BAO+hf(z). 

To sum up, after constraining the three models proposed in this work with the set of obser¬ 
vational data specihed in Sec. HV] it is seen hrst that, as a result of the joint effect of the 
vacuum and the particle creation rate, effective equations of state of phantom type can be 
obtained at present time and in the past. 


VI. COMPARISON WITH MODELS OF COUPLED DARK ENERGY 

At this point one may wonder whether models of DE featuring a nongravitational coupling 
with dark matter may also lead to an effective EoS of phantom type. As we will see, the 
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Figure 3. Left panel: la and 2a confidence contours for model II obtained from the joint analysis 
SNIa+GRB+BAO. Right panel: Evolution of the effective EoS in terms of the scale factor; the 
solid (black) line indicates the best-fit values, and the shaded region indicates the la uncertainty. 




Eigure 4. Same as Eig. [3] but using the set of data SNIa-|-GRB+BAO-|-Lf(z). 


answer is they very likely will not. 


The said models were proposed to lower the valne of the cosmological constant [33| and to 
solve or, at least, alleviate the coincidence problem (“why are the densities of dark energy 
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Figure 5. Left panel: la and 2a confidence contours for model III obtained from joint analysis 
SNIa+GRB+BAO. Right panel: Evolution of the effective EoS in terms of the scale factor; the 
solid (black) line indicates the best-ht values, and the shaded region indicates the la uncertainty. 



Eigure 6. Same as Eig. [5] but using the set of data SNIa+GRB+BA0+Lr(2;). 


and dark matter of the same order precisely today?”); see, e.g.. Ref. [3^. In any case, 
a conpling between the dark components seems natnral, though it must be small as it is 
severely constrained by observations -see 3^ for a review. 
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For a spatially flat, homogeneous, and isotropic universe, these models are characterized by 
the set of equations 


Pb + 3Hpb = 0 , 

(27) 

Pdm T Pdm Q ■} 

(28) 

Pde T 3Fr(l -|- W^Pde Q y 

(29) 

3H^ = 87lG{pb + Pdm + Pde) , 

(30) 


where Q stands for the energy transfer rate per unit of volume between the dark compo¬ 
nents and the EoS of DE is bounded by —1 < w < —1/3. This class of models implicitly 
assumes that the creation pressure either vanishes or is negligible. Because of the interac¬ 
tion, the energy density of DM, = nm, differs from the usual conservation expression, 
Pdm = Pdmo either because m varies with expansion or because the number density of 
DM particles does not obey the expression n = noa~^. 

Two distinct possibilities arise: (i) Q > 0, which means a continuous transfer of energy from 
DE to DM, and (ii) Q < 0 in this case, the transfer of energy would proceed in the opposite 
sense. For Q > 0, the effective EoS, tCg// = w + Q/{3Hpdm), results less negative than w, 
and thereby no phantom behavior can be reproduced (this holds irrespective of whether it 
is the number or mass of DM particles that varies). 

For Q < 0, one follows that iCe// can be less than —1. However this kind of model can 
be readily told apart from the ones discussed in the above sections. In particle creation 
models, DM particles are continuously added to this pressureless component, and thereby 
the amount of DM in the past is necessarily less, at any z > 0, than (at the same redshift) 
in the ACDM model that shares the same Dfeo and than the matter creation model. 
Likewise, in coupled models with Q < 0 the amount of DM in the past was certainly larger, 
at any z > 0, than (at the same redshift) in the ACDM model with identical Dm and ^Idmo- 
Consequently, the growth function, fg = dlnD+/dlna, in matter creation models is greater 
than in the ACDM and fg is also greater in the ACDM than in the said kind of coupled 
models. 

Further, as is well known, the transfer of energy from DM to DE worsens the coincidence 
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problem and violates the second law of thermodynamics 3^, and the conservation of qnan- 
tnm nnmbers conld be transgressed, especially if DE is the qnantnm vacnnm. 


In snmmary, models of conpled dark energy clearly differentiate observationally from creation 
models, and it is very nnlikely that they can reprodnce the EoS of phantom models. 


VII. CONCLUDING REMARKS 


Cnrrent observational data seem to favor an EoS for DE less than —1 (for a qnick snmmary 
of the state of art, see Fig. 1 in Ref. [^). If conhrmed by fntnre experiments, like Enclide 


37l |. we will confront the pnzzling sitnation that the cnrrent expansion stage of the Universe 


is likely driven by a phantom held, in spite of the fact that these helds are known to come 
with theoretical drawbacks of no easy hxing. 


In this paper we explored the possibility that the Wde determined by recent experiments is in 
reality an effective EoS that results from adding the negative EoS, Wc = —T/3H (associated 
to the particle production pressure from the gravitational held acting on the vacuum, j^. 


17|), to the EoS of the vacuum itself, wa = —1. Since the rate of particle production, E, 


is not known, we have assumed three diherent phenomenological Ansatze for the latter on 
grounds of simplicity, Eqs. ([H])-(IH]), and constrained their free parameters {(3 and ^dmo) of 
the corresponding cosmological models with the combinations SN la + GRBs + BAO and 
SN la + GRBs + BAO + II(z) of observational data sets. Here we wish to stress that we 
did not hx the value of the Hubble constant, Hqj at any point in the statistical analysis. As 
it turns out, models I and H suggest that the pressure of the vacuum combined with the 


particle crea 


Ion pressure helps explain that the EoS measured by recent experiments (see. 


e.g.. Refs. 2j- is less than —1. Model HI lies far away from explaining it. 


To sum up, the recently reported values of lower than —1 for the equation of state of dark 
energy may arise from the joint effect of the quantum vacuum and the process of particle 
production. This offers a viable alternative to the embarrassing possibility of helds that, 
among other things, violate the dominant energy condition, give rise to classical and quan¬ 
tum instabilities, and do not respect the second law of thermodynamics. 


17 









Obviously, phenomenological models of particle production different from the ones essayed 
here are also worth exploring. More important, however, is to determine the rate F using 
quantum held theory but, as said above, this does not seem feasible until the nature of DM 
particles is found. 
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